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The particle and event horizons are widely known and studied concepts, but the study of their properties, 
in particular their evolution, have only been done so far considering a single state equation in a deceler- 
ating universe. This paper is the first of two where we study this problem from a general point of view. 
Specifically, this paper is devoted to the study of the evolution of these cosmological horizons in an accel- 
erated universe with two state equations, cosmological constant and dust. We have obtained closed-form 
expressions for the horizons, which have allowed us to compute their velocities in terms of their respective 
recession velocities that generalize the previous results for one state equation only. With the equations of 
state considered, it is proved that both velocities remain always positive. 
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1 Introduction 

For any observer A we can define two regions in the instantaneous three-dimensional space t = to- The 
first one is the region defined by the events that have already been observed by A, and the second one is its 
complement in the three-dimensional space, i.e. the region that cannot be observed by A at a time to- The 
boundary between these two regions is the particle horizon at to, that defines the observable universe 
for A. Notice that the particle horizon takes into account only the past events with respect to A. Another 
horizon could be defined taking into account also the A's future. This horizon is the event horizon and 
it is defined as the hyper-surface in space-time which divides all events into two classes, those that will be 
observable by A, and those that are forever outside A's range of observation. This horizon determines a 
limit in the future observable universe [7]. 

A deep study of the horizons has been made before in [2] where only one state equation was considered. 
However, it is currently widely accepted that our present Universe is dominated by two state equations 
(cosmological constant and dust), one of them of negative pressure, that drives the universe into an accel- 
erated expansion. This fact leads us to make a deep study of more general situations that we summarize 
in a couple of papers, in this one, we make a deep study taking into account these two state equations, 
so that we can obtain results applicable to the currently accepted cosmology, and in the second one [5], a 
complete general study (at least from a mathematical point of view) has been made considering countably 
infinitely many state equations. 
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We are going to use the following cosmological parameters, for a flat concordance universe: 



Hubble constant 


Ho = 70.1 km s-^Mpc-i 


Current density of matter (dark and baryonic) parameter 


VLmo = 0.278 


Current density of cosmological constant parameter 


Q,Ko = 0.722 



This cosmological parameters have been obtained from [6] where a combination of data from WMAP, BAO 
and SNCONST is considered. Note that as a good first approximation, we are neglecting in this paper the 
parameter of density of radiation, f2rO, although it should be taken into account when considering times 
near the beginning of the Universe, when radiation dominates over dust, as indeed we did in [5]. 

This paper could be summarized as follows: in section one, we briefly introduce some basic cosmological 
concepts and establish the nomenclature to be used throughout the paper. In sections two and three, 
we define the formal concepts of particle and event horizons, and derive their integral expressions at 
any cosmological time. In section four we study the evolution of the horizons, obtaining more suitable 
expressions for the horizons at every time. Section five is devoted to gathering all the results obtained, 
discussion and conclusions. In this last section, some relevant graphics about the horizons and their 
derivatives in the at concordance universe are included. Throughout this paper, we use geometrical units 
where c = G = 1. 



2 Cosmological Prerequisites 
2.1 Proper Distance 

In order to introduce the particle and event horizons, we need to define the proper distance Dp, which is 
the distance between two simultaneous events at a cosmological time to measured by an inertial observer. 
Considering a homogenous and isotropic universe, we can write the Robertson- Walker metric, where a{t) 
is the scale factor (albeit with units of longitude in our system) and k the sign of the curvature: 

ds^ = -dt^ + a{tf (^ ^ ^^^^ dr'^ + r^dn'^^ (2.1) 

Isotropy of the Robertson- Walker universe allows us to consider any direction, i.e. 8, (j> constant and so 
dO.^ = 0. Therefore by the definition of the proper distance (where dt = 0) we have for a given time to: 



Dp^i ds^aoj^^^^ (2.2) 





where ao = a(to)- On the other hand, if we consider light rays {ds = 0), we have, from (2.1): 

dt _ dr p« dt' _ r dr _ Dp 

a ^ Vl - fcr2 ^ Jt a{t') ^ Vl - kr^ ^ 

Therefore the proper distance can be expressed as a distance measured using light as follows: 



(2.3) 



df' 

W) ^'-'^ 

So to speak, the previous formula represents the distance covered by the light between two points of the 
space-time, but considering an expanding universe since the a factor accounts for this expansion. 

2.2 Hubble Parameter 

The Hubble parameter H is defined as: 

ff=i^ (2.5) 
a dt ^ ' 

whose value at t = to is the Hubble constant Hq. The product of H times the proper distance, has 

dimensions of velocity and is known as the recession velocity: 

Vr = HDp (2.6) 
which physically is the instantaneous velocity of an object at a distance Dp with respect to an inertial 
observer. 

We are considering the universe as a perfect fluid with density p and pressure p, but it can also be 
approximated as composed by separate constituents of density pi and partial pressure pi , where all together 
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add up to p and p respectively. Furthermore, the i-th quantities are related by the linear state equation 
Pi = WiPi. Now using Einstein Field Equation for a perfect fluid, the Friedman equations, and considering 
the i-th state equation only, we can obtain [3, chap. 3] the i-th density in terms of the redshift z: 

= p,o(l + ^)'<'+™'' (2.7) 

Now we define the following magnitudes: 

3 2 

• Critical density as pc = —H 

OTT 

• Dimensionless energy density as = — 

Pc 

Pi 

• Dimensionless i-th energy density as ^li = — 

pc 

All this quantities, when referred to the current time, are denoted with a zero subindex. If we now substitute 
equation (2.7), pc and Qio in the definition of Qi, we obtain: 

n,{z) = (2.8) 

Finally, adding all the Qi (where the index i ranges over all possible state equations) leads to: 

H{z) = Ho 

3 Particle Horizon 



f^^(l + ^)3(i+».) (2.9) 

\ z = l 



3.1 At the present cosmic time 

As the age of the Universe and the light velocity have finite values, there exists a particle horizon Hp, that 
represents the longest distance from which we can retrieve information from the past, so it defines the past 
observable universe. Then the particle horizon Hp at the current moment is given by the proper distance 
Dp covered by the light coming from t — (the origin of the universe, where a hot big bang is assumed, 
compatible with the dominant equation of state) to to: 

Hp = \im Dp (3.1) 

t— ^0 

It is well known that the scale factor a is related with the redshift z by: 

l + z=^ (3.2) 

Combining equations (2.4), (2.5), the derivative of (3.2), and taking into account that z{to) = 0, it follows 
that the proper distance can be written as 

Particularizing equation (2.9) to a flat universe (il — 1) with non-relativistic matter (wm ~ 0) and cosmo- 
logical constant {wa = — 1) [3, chap. 3], we have: 

H{z) = HoV^M^ + ^y + ^Ao (3.4) 

Applying (3.1), (3.3), (3.4), and using the boundary condition that z — >■ oo when t — > 0, we obtain the 
expression of the particle horizon at the current cosmic time: 



dy 



(3.5) 



Ho Jo y/Qmo{l + y)^ + Qao 
3.2 At any cosmic time 

We have so far obtained the particle horizon for the current cosmic time (corresponding to 2 = 0). If we 
want to know the expressions of this horizon at any cosmic time, or equivalently at any other z' , we have to 
replace in (3.5) the constant parameters Ho, QmO and Qao (corresponding to z = 0) with the parameters 
corresponding to z' that we denote H{z'), Qrn{z') and Qa{z'). Where we recall from equation (2.8) that 
flrn{z') and ilA{z') are the dimensionless energy densities given by: 
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Q.a{z') =nAOjjJ^ 0.m{z) ^Q.mOjjJ^{l + z'f (3.6) 

Therefore, we obtain the particle horizon for any redshift z (where the prime is omitted in order to simphfy 
the notation): 

H{z) Jo ^nm{z){l + yr + nA{z) 



4 Event Horizon 

4.1 At the present cosmic time 

The event horizon represents the barrier between the future events that can be observed, and those that 
cannot. It sets up a hmit in the future observable universe, since in the future the observer will be able to 
obtain information only from events which happen inside their event horizon. According to its definition 
the event horizon can be expressed as: 

= lim (-Dp) (4.1) 

Proceeding as in section 3.1, but now taking into account that z — )■ —1 as t ^ oo, we obtain that the event 
horizon is: 

H ^ — r '^'^ (4 2) 

4.2 At any cosmic time 

Now proceeding as in section 3.2, the event horizon at any z is: 

He{z) = r '^^ (4.3) 

H{z) J_, ^r!„(z)(l + y)3 + nA(z) 



5 Study of the Evolution of the Horizons 

In this section we are going to obtain the expressions of particle and event horizons (as well as its derivatives) 
using hypergeometric functions. We would like to recall that not only is this approach interesting because 
we obtain analytical expressions for the horizons, but it is also a great advantage since these functions have 
been widely used and analysed. 



5.1 Mathematical Background 

The theory of hypergeometric functions can be found in many books devoted to advanced calculus. For 
general theory see for example [1] or [4, chap. 9], the latter is the one we are following here. In what follows, 
we are going to establish the minimum amount of theory that we need for our purposes. 



Definition 5.1. 

2-Fi has the following (integral) definition: 



' r(/?)r(7-/3)yo i^-tx)- \arg{l-x)\<n ^^'^^ 

where T{z) is the gamma function, defined by: 



r(z) = / e-'f'-'dt, Re{z) > (5.2) 
Jo 

and by analytic continuation for every z G C apart from 2 = 0, —1, —2 . . . where simple poles appears. 



Remarks. 

5.2 Let us recall the important factorial property of the Gamma function, V[z + 1) = zV{z) for every 
z G C apart from z = 0, —1, —2 . . . and r(l) = 1. 

5.3 Since no distinction between hypergeometric functions pFq is necessary, from now on, we are going 
to write 2F1 = F. 
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5.4 The function F, as a function of x, is defined in the whole complex plane cut along [l,oo), though 
we would only need real negative values. 

5.5 The function F can be extended for every a, j3 £ C and for every 7 7^ 0, —1, —2 . . . through recursive 
equations, in fact, F{a, P;"/; x)/r{'y) is an entire function of the parameters a,/3,7. For more details 
see [4, sect. 9.4]. 

5.6 When we restrict F, considered as a function of x with fixed real parameters, to real values of i, we 
obtain a real- valued function F : (— 00, 1) — >• R. 

5.7 Applying the properties of the beta function B{x, y) 

B{x,y)= f f"-'(l-t)^-^di, Re{x),Re{y) >0 

r(7) _ 1 







we can rewrite 



r(/3)r(7-/3) B(/3,7-/3)- 



Now we are going to state some useful properties (where we are always considering that the third argument 
is different from 0, —1, —2 . . . and the fourth one verifies \arg{l — x)\ < n). All of them can be found in [1] 
or in [4, sect. 9.2,9.3,9.8,9.5]. 

Properties. 

5.8 F(a,/?;7;0) = 1 



dF ■7—1 
5.9 ^(a,/3;7;:c) = -^^ 
ox X 



F{a,p-r,x)-F[a,p--i-l,x) 



5.10 hm F(a,/3;7;a;) = ^ Re{^ - p - a) > Q 

5.11 F(a,/3;/3;s) = (1-a;)-" 



5.12 F{a,p-r,x) = ^ . F(q,7-/3;7 



l~xY V ' x-\ 

5.13 hm na,A7;-a;) =0 ""'^L^ 

The last expression follows by one of the analytic extensions of F given in ]4, sect. 9.5] together with 
property 5.8. Now we are going to obtain some particular values of F that we will need. 



5.14 / ^ 2F\\}--\\-A\ if A>0 



Lemmas. 

f- 00 

/i x/i^TA V2'6'6' 

da 1 „ fl 1 4 1 



5.15 / = —F ) if A>0 

Jo Vs^ + A V3 V2 3 3 



5.16 -j;-{a,p-p + l-x) = -^ 
ox X 



F{a,p-p + l-x)- 



(1 - xY 



5.17 lim F ( ^, 1; 1 + c;a; ) = 00 if 0<c<| 
x-^i^ \2 J 

Proof. 

1. Beginning from the right hand side of the equation, substituting in the definition (5.1), applying the 
factorial property of the Gamma function and then making the change of variable s = t~^^^ leads to 
the left hand side of the equation. 

2. Analogous to the previous case, but now making the change of variable s — t^^^ . 

3. This statement follows immediately from properties 5.9 and 5.11. 
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4. Notice first that property 5.10 cannot be applied as the conditions on the parameters do not hold. 
The result follows applying property 5.8 to the analytical extension in |2; — 1| < 1 cut along [l,oo) 
that can be found in [4], equations (9.5.7) or (9.5.10). The hypothesis of this equations are satisfied 
as c G (O, I) and a; € K. ■ 



It is important to note that in the third equation, we have managed to established the derivative of F (with 
one restriction in the parameters) in terms of the F itself. Finally, we introduce the auxiliary functions 
A, B and gather some straightforward computations that we will need to obtain the derivative of both 
horizons. 



1 



B ; 



»a(^) 

^m.{z) Q,mQ (1 + z)'-^ 
1 «™0(^ + ^^3 



A Qao 
H{z) = ffo\/a^\/BTT 
dz ao da ao 



dA 
dz 

dB 

dz 



3A 

1 + 2 

jiB_ 

1 + 2 



dt 



a2 dt 



H = -{l + z)H{z) 



(5.3) 

(5.4) 

(5.5) 
(5.6) 



Where the last expression is obtained applying first equation (3.2), then equation (2.5) and finally equation 
(3.2) again. 



5.2 Expressing the Particle Horizon through Hypergeometric Functions 

The expression of the particle horizon (3.7) can be expressed as follows, where z should verify z £ (—1, oo) 
according to equation (3.2): 



H,{z) 



dy 



H{z)Jo Vn™(2)(l+J/)3 + nA(2) H{z)y^aJzjJl s'^ + A{z) 



ds 



I f[^- i.I.-A(2)^ 



? f(1 l-l-^Mz^ 



Finally, by definition of A (5.3): 



HAz)= ^^Fri,l;^;-A(2) 



(5.7) 



5.3 Obtaining the derivative of the Particle Horizon 

Now we are going to obtain the derivative of Hp. Notice that in the expression of Hp, the parameters of F 
verify 7 = /3 + 1, so we can use lemma 5.16 and obtain F with the same parameters. For that reason and 
in order to simplify the notation, we omit the parameters and we also omit the argument 2 in the A and 
B function. 



dHp _ 1 

dz //oV^mO 
1 



Ha^/0.mO 
1 



dA dz , r-, dA\ dF , 

2^^F i-A) + 2VA — {~A) 

2^/A \ dz J dx ^ 



3^Fi-A) + ^jFi-A)- 



1 + 2 
2VA 



(1 + 2) HoVu:::^, 



Fi-A) 



1 + 2 

1 



ViTA 



B^l/A 



eg. (5. 7), (5. 5) 



(1 + 2) HoV^iVTTB 



^ Hpjz) _ 1 

1 + 2 {l + z)H(z) 

Finally, applying the chain rule and equation (5.6), we obtain: 



^=Hp{z)H{z) 



(5.; 
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Note that Hp{z)H{z) represents the recession velocity of the particle horizon (2.6). Its physical meaning 
is explained in subsection 6.5, where all the conclusions are provided. 



5.4 Expressing the Event Horizon through Hypergeometric Functions 

Analogously for He we have: 



He{z) 



dy 



a — l+y 



ds 



H{z) ^n„(2)(i + yf + n^iz) H{z)y/a:4zj Jo V^^^Ta 

1 



1 



1 4 1 



F 



114 1 



H{z)^n;4^)VA V2'3'3' AJ H{z)^nI(Y) V2'3'3' A 

Where in the last equality, the definition of A (5.3) is used. Now using equation (3.6) we obtain the 
following expression: 



He{z) = 



(5.9) 



5.5 Obtaining the Derivative of the Event Horizon 

We are again in the hypothesis of lemma 5.16, and then it is meaningful to omit the parameters. 
dH, 1 f dB\ dF 



dz HoV^Ao V dz J dx 
1 1 



(1 + z) HoV^ 
HJz) 1 



F{-B) + 



eg. (5. 9), (5. 5) 



1 + z HoV^VlTB 



1 + 2 {l + z)H{z) 
Finally, applying the chain rule and equation (5.6), we obtain: 



dHe 
dt 



= He{z)H{z) - 1 



(5.10) 



Where now He{z)H{z) is the recession velocity of the event horizon. Its physical implications will be 
explained in subsection 6.5. 



6 Results and Interpretation 
6.1 Some relevant values 

Below we show some important values of the horizons, summarized in table 1, and the required computa- 
tions to obtain them. 





z 


t 








He{z) 


dHp 

dt 


dH, 
dt 


Origin of the universe 


oo 





1 











3c 


oo 


Current time 





to 


0.278 


0.722 


2.672 c 


0.958 c 


4.409 c 


0.128c 






Future 


-1 


oo 





1 


oo 


C 

i/oV^AO 


oo 






Table 1 - Some important values of the horizons and their velocities in physical units. Notice that 
the numerical values appearing on the current time values, depend on both Qao and QmO- 
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Hp{oo) - Notice that A — >■ when z — >■ cxj and use property 5.8 in (5.7). 
He(oo) - Notice that B — )■ oo when z — >■ oo and use property 5.13 in (5.9). 

^^^(oo) - Again, A — >■ when z oo, then apply property 5.12 to eq. (5.8), which leads to: 
at 

(oo) - Again, B — > oo when z oo, then apply property 5.12 to eq. (5.10), which leads to: 

dHe r.f^ . . 1 B \ , 5.17 

= F[ -, 1: 1 + -; — - 1 > oo 

dt V2'' 3'-B + iy 

Hp{-~1) - Notice that A ^ oo when z — )■ —1 and apply property 5.12 to eq. (5.7), which leads to: 



2 A ^[1 ^ ^ 1 A \ 5.17 
Hv = ; a/- F -, 1; 1 + - ; > oo 

• He{~l) - Notice that B — > when z — > —1 and use property 5.8 in (5.9). 

• To determine the values of the derivatives for z = —1, just substitute H{—1) = Hoy/ilAo in equations 
(5.8) and (5.10), and use the values Hp{—1) and He{ — 1) we have just obtained above. 

• Finally, we obtain the values at the current cosmic time through numerical calculus using (5.7), (5.9), 
(5.8) and (5.10). 

6.2 The Velocities are always Positive 

It is clear from equation (5.8) that the velocity of the particle horizon is always positive, but from equation 
(5.10) we cannot conclude the same for the event horizon. Let us see that with a little more effort we 
are able to obtain some information about the behaviour of the horizon velocities, including the fact that 
they are indeed always positive. In order to do that, we compute the second derivative which is quite 
straightforward using all the previous computations: 



df2 \2{1 + B) P ' y dt^ \2{1 + B) 

If we equal both equations to zero, we obtain (where the correspondent parameters have been omitted): 



Fori?,: F(-A).^ Fori/.: Fi-B) ^ '-^^ 

From the definition of F, it is easy to prove applying elementary inequalities rules, that when x £ (— oo, 1) 
with positive parameters the function F{x) is strictly positive and increasing. Then on one hand, we have 
that F{—x) is positive but decreasing when x £ (—1, oo).^ Computing the derivative with respect to A of 
the right hand side of first equation, we obtain on the other hand, that it is increasing where it is defined, 
takes positive values in the interval / = (—1, 1/2) and negative for larger values. Gathering these two facts 
we have that if there exists a solution Ao for the first equation, it must be unique and Ao £ i. Clearly 
when A — 0, the right hand side is 1 and from property 5.8 i^(0) = 1, so the unique solution is Ao = 0. 
Analogously for the second equation (but now I = (—1,2) that does not affect the argument) we obtain 
that i3o = is the unique solution for the second equation. 

Therefore, from the definition of A and B, we have that the acceleration of Hp is zero only at the origin 
of the universe 2 — > oo and the acceleration of He is zero only at the far future 2 = — 1, and for the rest of 
the values of 2, both accelerations are positive. 

As at the origin of the universe both velocities were positive, we conclude they are always nonnegative. In 
fact, they can only vanish for the limits 2 — > 00 and 2 — >■ — 1, and according to table 1, it only happens at 
the far future for the event horizon speed. We have included on the appendix (page 11) some figures that 
show the behaviour of both horizons and their derivative. 

and B, that by definition are both positive when z ranges over (—1, 00), would play the role of the x. 
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6.3 Limit to just one state equation 

In this section we are going to prove that in fact we are generahzing the previous results obtained taking 
into account just one state equation such as the ones stated in [2]. In general, it can be proved that for a 
unique state equation the following results are obtained for any cosmic time (e.g. using (3.1), (3.3), (2.9) 
and (2.8) for just one state equation, and (5.6) to compute their derivatives): 

\ 2 dH; 3{l + w) 

^-(") = ]?^MTT3^ ^ = TT3^ if«;>-i/3 

Where the * stand for the fact that we are taking into account just one state equation, in contrast with 
the equations obtained in this paper taking into account two state equations. 

• If it dominates the state equation of matter, we have w = Q, then there only exists the particle horizon 
whose equations, as proved in [2], are: 

. , . 2 dH* 



In this case Q,m{2) = fimo = 1 and Q,k{z) = Qao = 0, and therefore H*{z) = Ho\/ (1 + 2)^. Now 
substituting the A function in equation (5.7), we have (omitting once again the parameters): 



HAz)- „ z, F{-^j-^]^-- ^ - F{0) 



5.8 



i/0v/fim0(l + 2)3 V f^m0(l + 2)V Ho^/(TTW ^* W 

'-§=Hi.)HA.) + l^H*i.)^^+1^3 

If it dominates the state equation of cosmological constant (i.e. w = ~1) then there only exists the 
event horizon whose equations are: 

_ _J_ ^ _n 

H*{z) dt 

In this case rimiz) = flmo = and ^a{z) = ^ao = 1, and therefore H*{£) — Hq. Then, equation 
(5.9) leads to: 

1 ^/ S^mO /I , x3\ 1 rn 5 S 1 1 



H.iz)= - F -^(1 + ,)^^^F(0) 



HoV^Uo V ^Ao J Ho ' ' Ho H*{z) 

=//(.)//.(.)- 1 = H*(.)^ -1 = 
So we conclude that indeed, our results generalize the previous ones obtained in [2]. 

6.4 The Case of the Cosmic Background Radiation 

The cosmic microwave background radiation that we measure today comes from a spherical surface called 
the surface of last scattering. In fact, it is not a surface but it has some thickness. In the current 
Standard Model, the surface of last scattering is at redshift zis ~ 1089 with a thickness of Az — 195 [8]. 
The proper distance to the mean value of this surface is: 

1 r^'^^^ dv 
Dp = — / ^ ^ = 14 086.37 Mpc 

Ho Jo ^ri„o(l + y)^ + f^AO 
The values of particle and event horizons at the current cosmic time are according to table 1 (with the 
numerical constant we mentioned at the beginning of the paper): 

Hpo = 14 577.72 Mpc Heo = 4 823.66 Mpc 

We notice that from the integral definition of the particle horizon (3.5), the distance to the surface of last 
scattering is always smaller than the particle horizon, which agrees with the experimental data. The light 
which is reaching us at the current moment from objects located over the surface of last scattering was 
emitted at a time corresponding to a redshift of zis- As time goes by, the distance to the surface will be 
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higher, but always less than the particle horizon, so in the future we get information emitted from the 
surface of last scattering. 

On the other hand the distance to the last scattering surface is greater than the event horizon, which means 
that the light emitted at the present moment from this surface, is never going to reach us. It does not 
mean that in the future we are not going to see background radiation, but that the radiation we can see 
at the current time and which we will see in the future is the one inside the present event horizon. This 
can be generalized to any object in the universe, and drives to the conclusion that in an ever accelerated 
universe (i.e. dominated by a negative pressure) with several state equations, where the initial expansion 
phases were dominated by positive pressure terms, although the particle horizon guarantees that an object 
can be seen, the event horizon prevents to see photons beyond a time given by the instantaneous event 
horizon at that time. Summarizing, the particle horizon defines the events that could be observed at a 
given cosmic time, while the event horizon defines the events that will be observed in the future. 

6.5 Conclusions. Recession of the Cosmological Horizons 

In this paper, we have first obtained analytical expressions for the particle and the event horizons consider- 
ing two state equations (cosmological constant and dust as accepted in the concordance cosmology). They 
have allowed us to compute their first derivative, obtaining these extremely simple equations: 



where HpH and H^H are the recession velocities of the particle and event horizons respectively. These 
expressions generalize the previous stated results [2] derived for just one state equation. Notice that all 
these results and values are independent of the current values of the involved constants. 

The equations of the velocity of the horizons have a remarkable physical meaning. As the recession velocity 
HpH is the instantaneous velocity of an object located at the distance of the particle horizon Hp, from 
the first equation we deduce that the instantaneous velocity of the surface of the horizon particle is faster 
by the speed of light c (in our units c = 1), than the one of the objects over the particle horizon, then 
more and more objects are entering into the particle horizon and they will never get out. Analogously, the 
second equation stands that the recession speed of the event horizon is slower by c than that of the objects 
over the event horizon, and then more and more objects are disappearing from the event horizon and they 
will never get into again. 

Finally, computing the second derivatives of both horizons we have proved that under our assumptions, 
their velocities remain positive throughout the history of the universe. 
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dHp 
dt 



Hp{z)H{z) + 1 



dH, 
dt 



H,{z)H{z) - 1 
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A Graphical Behaviour of the Horizons 



In this appendix we include the graphics of the horizons and their velocities obtained with numerical 
calculus. As the hypergeometric functions are really well implemented in almost all mathematical programs, 
obtaining these graphics is extremely easy. 

A.l Particle Horizon 
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Figure 1 - Particle Horizon. At the origin of the Universe (right hand side of the graphic on the 
right), the particle horizon tends to zero, then it increases rapidly with time tending to infinite at the 
far future (left hand side of the other graphic on the left). The dot in both graphics represents the 
connecting point between them and corresponds to the current time 2 = 0. 
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Figure 2 - Velocity of the Particle Horizon. At the origin of the Universe it tends to the constant 
value 3c ~ 2.91 ■ 10""Mpc/s and increases rapidly to infinite, just as the particle horizon itself. 
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A. 2 Event Horizon 



5050 



5000 



4950 



4900 



4850 



4800 



Future Event Horizon 



Past Event Horizon 




5000 



4000 



3000 



2000 



1000 



-0.6 -0.4 

z - Redshift 




■ Redshift 



Figure 3 - Event Horizon. At the origin of tlie Universe, the event horizon tends to zero, then it 
increases with time, tending at the far future to the constant value c (//o V^^ao) — 5 033.08 Mpc 
with no velocity i.e. with zero slope. 
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Figure 4 - Velocity of the Event Horizon. At the origin of the Universe it tends to infinite and 
decrease to zero at the far future with no acceleration (again with zero slope). 
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